Colorful SU(2) center vortices in the continuum and on the lattice 

Thomas Schweigler, Roman H611wicsci00 and Manfried Fabei@ 
Institute of Atomic and Subatomic Physics, Vienna University of Technology, Wiedner Hauptstr. 8-10, 1040 Vienna, Austria 

Urs M. Heller 

American Physical Society, One Research Road, Ridge, New York 11961, USA 
(Dated: February 12, 2013) 

The spherical vortex as introduced in [Phys. Rev. D77, 014515 (2008)] is generalized. A continuum 
form of the spherical vortex is derived and investigated in detail. The discrepancy between the 
gluonic lattice topological charge and the index of the lattice Dirac operator described in previous 
papers is identified as a discretization effect. The importance of the investigations for Monte Carlo 
configurations is discussed. 

PACS numbers: 11.15.Ha 



a Funded by Austrian Science Fund (FWF) under contract P22270-N16. 
t hroman@kph.tuwien.ac.at 
* fabcr@kph.tuwicn.ac.at 



2 



I. INTRODUCTION 



The distribution of topological charge density is closely linked to the phenomena of the axial anomaly and sponta- 
neous chiral symmetry breaking. There are stron g hi nts that the QCD vacuum is dominated by center vortices ■ 
Center vortices can explain color confinement [7H10I] and seem to be also of paramount importance for spontaneous 
chiral symmetry breakin g [ll| - [l5| . Center vortices can contribute to the topological charge density through intersec- 



tion and writhing points ll|, but also through their color structure. The prototype of the later contribution in SU(2) 
gauge theory is the spherical vortex introduced in previous articles of our group |1 6l — ll 8I | . In this paper, a generalization 
of the original spherical vortex is constructed. Subsequently, the generalized spherical vortex is investigated in detail 
on the lattice and in the continuum. 

In section HH we start with a description of the original spherical vortex on the lattice as introduced in [l(| . From 
this lattice gauge configuration, we then derive the corresponding continuum object. With the continuum form, a 
more general spherical vortex can be described. The action and topological charge density of such generalized spherical 
vortices in the continuum are discussed in section IIIII and IIV1 We show that the large action of the original spherical 
vortex encountered in previous investigations can be significantly reduced by spreading the vortex over several time 
slices. In section [Vj a generalized spherical vortex on the lattice is derived from the continuum form. Subsequently 
its gluonic and fermionic properties are investigated. In particular, we show that the discrepancy between the gluonic 
lattice topological charge and the index of the lattice Dirac operator seen for the original spherical vortex is simply a 
discretization effect. We conclude with remarks on the contributions of colorful center vortices to topological charge 
and chiral symmetry breaking in Monte Carlo configurations. 

II. THE ORIGINAL SPHERICAL VORTEX 

The original spherical vortex on the lattice as introduced in [f(| is given by the lattice links 

Ui(x) = 1 

cos [a (\f — ro\)] 1 — i e r ■ a sin [a (\f — ro\)] for t=l (1) 



1 else 



Here a; = {:pi, X2, £3, t} stands for the lattice site and r denotes the spatial components only. The unit vector pointing 
from the vortex midpoint r*o to the spatial lattice site r is denoted by 



-. _ r -r 
\r — fo\ 

Ufj,(x) is the link connecting the lattice site x with the site x + /t, where p, is the unit vector pointing in the \x 
direction. The gauge links on the lattice represent the adjoint gauge transporters and, therefore, transform like 
Up{x) — > fl(x + fi) J7 M (x) fl(xY under a gauge transformation Q(x). Indices denoted by Latin letters run from I to 3 
(spatial indices), while indices denoted by Greek letters run from I to 4 (spatial plus temporal indices). The function 
a(r) rises linearly from — ir to 0, i.e., 

' 1 for r < R - I 

{^±p: {orR-l<r<R+%. (2) 
for r>R+ \ 

Note that the spherical vortex configuration posses only non-trivial temporal links U4. All spatial links Ui are 
trivial and therefore the gluonic lattice topological charge Q = 0. On the other hand the index of the Dirac operator 
ind [D] = ri- — rt+, defined as the difference between the number of left- and right-handed chiral zeromodes, is always 
I. This gives a discrepancy between the index and the gluonic topological charge. Note that no violation of the 
index theorem ind [D] = Q occurs as the theorem is only stated for the continuum limit. When cooling, the gluonic 
topological charge approaches 1 and the discrepancy ind [D] ^ Q gets resolved fill - fisj . 

Naively, one would perform the continuum limit of the spherical vortex by simply splitting the temporal links and 
leaving the spatial links trivial. This would give gauge configurations of the form 



Ut(x) = 1, 



Ui(x) = cos [w(t) a (|r — ro|)] 1 — i e r ■ a sin [w(t) a (\r — ro|)] with ~^^w(t) 

t 



(3) 
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where w{t) determines how the temporal links are split. As one can easily see, a singularity occurs at the vortex 
midpoint when ^ w(t) ^ 1. Moreover, the gluonic topological charge Q clearly remains independent of w(t) which 
doesn't agree with the cooling history of the original vertex, Eq. (pj. Therefore, this simple procedure seems not to 
be suitable for performing the continuum limit of the lattice spherical vortex. 

In order to find a better continuum limit and to gain understanding of the spherical vortex, one has to apply a 
gauge transformation to the link configuration given in Eq. ([1} . With the lattice gauge transformation 

for Kt<t 
W [1 else 

where 

g{f) = cos [a (|r — ro|)] 1 + i e r ■ a sin [a (|r — ro\)] , (4) 

the spherical vortex becomes 



Ui(x) 
U 4 {x) 



g (r + ij g (r) t for 1 < t < t g 



g (r)^ for t = t g 
1 else 



else / 

(5) 



From this it becomes clear that the spherical vortex represents a transition (in the temporal direction) between two 
pure gauge fields. The transition occurs between t = 1 and t = 2. No transition occurs between t = t g and t = t g + 1 
since all plaqucttes there are trivial. For t < 1 the pure gauge field is trivial, i.e., the winding number n w — 0. For 
t > 1 the pure gauge field is generated by the hedgehog gauge transformation, Eq. @, which has n w = 1 [2(|. This 
means that in the continuum limit, the spherical vortex (representing a transition between a vacuum with winding 
number and 1) has topological charge Q = 1 (20| . It can be regarded as a squeezed instanton with center vortex 
structure. 

Assuming an infinitely big temporal extent of the lattice and taking t g — >■ oo, the continuum field corresponding to 
Eq. ^ can be written as 

A ll = i(t) i(dtf)gi, (6) 

where g is given in Eq. Q and f(t), determining the transition in temporal direction t, changes from to 1 within 
one lattice unit. Clearly, one could also use a function f(i) that changes more slowly. In this case, the corresponding 
lattice object cannot be represented by configurations of the form of Eq. ([1]) anymore. We refer to the construction 
with a smoother function f(t) as a generalized spherical vortex. One construction of such a generalized spherical 
vortex on the lattice will be investigated in section [Vj Before doing so, we discuss the action and topological charge 
density of the continuum gauge field, given in Eq. (JS]) in the next two sections. 

III. ACTION OF THE SPHERICAL VORTEX IN THE CONTINUUM 

In this section, we calculate the action for the generalized spherical vortex in the continuum given in Eq. ((5]). For 
simplicity, we set the vortex midpoint ro to in this and the next section. Moreover, we use the notation |r| = ?\ 

To explicitly evaluate the action S for the given in Eq. we have to evaluate i (d^g) g^ where g is given in 
Eq. (Q| . To do so, we first have a look at i (d^g) g^ for a general gauge transformation 

g = qa (JQ+iq - a, 

where go and q are constraint by q$ + \q\ 2 — 1. Using the identity <j a a b = S ab 1 + ie abc a c and 

qad^qo + qkd^qk = ^ (<7o<7o + 9fe9fc) = \d^ L (1) = , 
one gets, after a few lines of calculation, 

* (dp.g) 5 f = s • [(d^qo)q- qo (d^q) + qx (d^q)} . (7) 
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For the g given in Eq. (j4|) we have go = cos(a(r)) and q '= e r sin(a(r)). Inserting this into Eq. ([7]) and multiplying 
withf (t) gives (A„ = 

A? = f (t) ■ 2 • ^^p- sin (a(r)) cos (a(r)) + ei a /c £fc sin 2 (a(r)) 

1 ZiX a \ ( 8 ) 
- (5 ia - sin (a(r)) cos (a(r)) ^""'M 

^ = 0. 



In the following, we use the notation 

A- = f(i)A£ + . (9) 

The field strength tensor is given by 

= M - - e abc AlAl . 

The different elements of the field strength tensor can be simplified considerably for our particular gauge field config- 
uration. Let us first have a look at the elements of the form Ffa. Using the fact that the gauge field has no temporal 
components, i.e., A\ = 0, we can simplify these elements to 

K = d 4 A? . (10) 

The Fjj can be simplified too. With the notation ([9]) we get 

F% - f (t) diA a j+ - f (t) - f (t) 2 e abc A b + A c j+ 

= f (t) - d 3 A« + - e abc A b + A^ + (f (t) f (t) 2 ) e a6c 4 + ^ c + (11) 

= f(t)(l-f(f)) t ahc A\ + A] + . 
Using these results, we can now evaluate the action. Let us first evaluate trc [F^F^]. With (fit))) and ((SJ) we get 

1 ™ ™ 1 / d 



tr c [^4^4,] = ^ *X*X = 2 ( ) A i+ A i+ 

Inserting the explicit form of Af + one can evaluate 

^■ + =4^sin(a(r)) 2 +a'(r) 2 

Similarly we get 

tr c [FijFiA = \ F%F« = f (i) 2 (1 - f (t)) 2 e abc A b + A c j+ e ade Af + A e j+ , 
where Eq. (fTTj) has been used in the last step. With the help of a computer algebra program, one can evaluate 

e abc A b + A c j+ e ade Af + A e j+ = psin(a(r)) 2 (sin (a(rjf + 2rV(r) 2 ) . 
Combining these identities and switching to polar coordinates gives 

= 2^ ( / dt i (1 - f (t))2 / dr sin ( a W) 2 ( sin ( Q W) 2 + 2r 2 a'(r) 2 ) 



(12) 



+ J dt(j t i(t)\ J dr 16tt (2 sin (a(r)) 2 + rV(r) 2 ) 
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So far, the calculation has been done for a general f (t). Now we will choose f (t) as the piecewise linear function 

!0 for t < 1 
for 1 < t < 1 + At , (13) 
1 for t > 1 + At 

where At stands for the duration of the transition. Inserting f At(i) into Eq. (JT2J) gives 




One can see from Eq. (|14p. that as long as none of the spatial integrations gives zero, the action diverges for both 
At — > oo and At —> 0. Note that the first term in Eq. ([LI"]) represents the magnetic and the second term the electric 
contributions to the action. This means that for At —> 0, the action is purely electric, and for At — > oo, it is purely 
magnetic. In between, there is a minimum where the electric and magnetic contributions are equal. 

Performing the spatial integration is a difficult task, which was done with a computer algebra program. The result 
for a general ratio of d/R is lengthy and therefore not explicitly given here. However, let us state some properties of 
the result. If we minimize the action for a given R and d with respect to At, we get an expression that depends only 
on the ratio d/R. This expression falls monotonically with d/R in the allowed range < d/R < 2. It diverges for 
d/R — > and approaches its minimum for d/R 2. The absolute minimum of the action (minimized with respect 
to At, R and d) is given by S m i n = 1.667 Sinst {R = d/2 « 0.305Ai). Here we denoted the action of one instanton 
by Sinst = l/(2<? 2 ) 167r 2 . This action serves as a lower bound for objects with topological charge \Q\ = 1. Fixing the 
ratio d/R gives the action as a function of R and At. For d/R = 1 we get 

S(At) _ 0.4358 Af + 3.722 R 
Sinst R At 



IV. TOPOLOGICAL CHARGE DENSITY OF THE SPHERICAL VORTEX IN THE CONTINUUM 

In this section, we calculate the topological charge density q(x) of field configurations of the form of Eq. ((6]). For 
now, the calculation is done for a general gauge transformation g and a general f (t) with the restriction, that they 
have to be chosen in such a way, that the resulting fields are free of singularities. Moreover, we will assume that g is 
independent of the temporal coordinate t. In the following, we again use the shorthand notation Ai+ = i (dig) g' . 

As is well known [2(|, q(x) can be calculated as the full derivative 

q(x) - O^K^x) , with K a = (KdpAZ, - ^e abc A a a A b A c y 

Let us now split d^K^ into a spatial (diKi) and a temporal part (c^-fQ). First, we treat the spatial part. With 
the assumption that g is independent of the temporal coordinate, one can easily see that the gauge field © has no 
temporal component, i.e. Ai — 0. Therefore, Ki evaluates to zero: 

/ . n ^ m 1 abc a a Ab * ^ \ 



16tt 2 V 7 3 J 16tt 

e ijk A^Al = - ^3 e ijk A a j+ A% + f(t) f'(t) = . 



16tt 2 16tt 2 
This means, that we can identify the topological charge density with the temporal derivative of K^, i.e., 

q(x) = d^Ki . 

From the definition of K u given above, we see that we can write K4 as 



K 4 = - 

16 



l — U{tf e ijk Al + djA% + - f(t) 3 l - e ijk e abc A% + A b j+ A k _^ , 
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where we have used eujk — — £ijk- Using the fact that the field strength vanishes for pure gauge fields, i.e., 

Fly = d*K+ - d » A l+ - eabcA l+ A i+ = °> one can rewrite the term e Hk A t+9jA% + as | e ijk e abc A a i+ A b j+ A% + . Using 
this result, we get 



K ^ = ~\Mi\\2 m 3 m ) Cijk 6 A "+ A j+ A k + ) • (16) 

The topological charge density is simply the temporal derivative of this expression. Note that ey^ e abc Af + A b + A^. + 
is proportional to the winding number density of the gauge transformation. This means that the spatial dependence 
of the topological charge density is given by the winding number density and the temporal dependence is given by 
(1 - f(f)) f(f) f'(t). With the g defined in Eq. gj), the a{r) defined in Eq. © and the f(t) defined in Eq. the 
topological charge density evaluates to 

3 , Mr-fl) \ / 1 _ t 2 \ fl forf<t<I + At and i? - f < r < i? + f 
91 ' j ~ dTrr 2 ° S 1, d / \l~Al At 3 ) \0 else ' ( ' 

As can easily be checked, integrating this expression (r 2 dr and dt) yields Q = 1. 



V. THE GENERALIZED SPHERICAL VORTEX ON THE LATTICE 



We now put the generalized continuum spherical vortex of Eq. ^ onto the lattice with periodic boundary conditions. 
Clearly, the field as given in Eq. (JSJ does not fulfill periodic boundary conditions in the temporal direction. Af, vanishes 
for t — > -co (and also for r — > co) but not for t —> oo. If we want to get a vanishing gauge field A^ for t —> oo by 
gauge transforming the field of Eq. we have to find a gauge transformation that equals 1 for t — > — oo and g* for 
t — y co. As can be shown easily by continuity arguments, there is no continuous gauge transformation that fulfills 
that criteria for the g given in Eq. (H|). However, there is still the possibility to put the field onto a lattice of infinite 
size. Then, one can transform the links for t — > oo to unity and subsequently close the lattice by periodic boundary 
conditions. Performing this procedure for the continuum field of Eq. ([6]) with f(i) given in Eq. (|13p yields 



Ui(x) 



U 4 {x) 

The functions g{f) and a(r) are again given by Eqs. (U]) and ©. 

Let us now have a look at the gauge action and topological charge of this generalized spherical vortex on the lattice. 
The gauge action as function of the temporal extent At is plotted in Fig. [TJ,) . Note, that the action of the lattice 
object matches the action of the underlying continuum object pretty well. However, it systematically underestimates 
the continuum value. In Fig. [TJj) we show the topological charge density as function of At for three values of R = d. 
As one can see, the discrepancy between the topological charge Q on the lattice and in the continuum (which is 
always I) is dramatic for small temporal extent At of the spherical vortex. For At = 1 the topological charge of the 
continuum object is not recognized. However, one still gets a zeromode for the corresponding lattice Dirac operator, 
i.e., the fermions still see the topological charge of the underlying continuum object. This is not unexpected as the 
vacuum to vacuum transition is still present in the lattice object. To be more precise, the lattice samples the field 
before and after the transition. The transition itself falls between two time-slices. As discussed in section IIV1 the 
transition carries the topological charge. Therefore, by missing the transition, one also misses the topological charge. 
From Fig. Q] we see, that also for very big At (slow vacuum to vacuum transitions), the lattice topological charge 
doesn't quite approach one. However, increasing R, keeping R = d in the example shown, increases Q towards one. 
Since R is the only spacial scale, 1/R acts as the spacial lattice spacing, showing that the deviations from Q = 1 at 
large At (At > 10) are discretization effects in the spatial directions. 

To conclude this section we briefly discuss how the fermionic results depend on the temporal extent of the vortex. 
For the fcrmion fields we use antipcriodic boundary conditions in the temporal direction and periodic boundary 
conditions in the spatial directions. With these boundary conditions, the Dirac operator possesses always exactly 
one zeromode for gauge configurations of the form of Eq. (|I8[) . independent of the values of the parameters. This 




(t-l)/Ai 



for 1< t < 1 + At 
for I + At < t < t g 
else 



(18) 



g(r)^ for t = t g 
1 else 
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FIG. 1. a) Comparison of the gauge action (in units of the instanton action Sinat) on the lattice, (blue) diamonds, with 
the gauge action of the continuum object, (red) line, given by Eq. (|15[) as function of the temporal extent At of the spherical 
vortex. The radius as well as the thickness of the vortex are given by R — d = 6. The midpoint of the vortex is half integer, 
b) Comparison of the gluonic topological charge Q(At) for three values of R = d. The lattice results (calculated with the 
plaquette definition) are shown as (green) squares, (blue) diamonds and (black) circles. In the continuum, Q(At) = 1 for all 
At as shown by the (red) horizontal line. The lattice sizes N% x N t are chosen to fit the spacial and temporal extents of the 
spherical vortex. 



zeromode is more or less located at the vortex. For the calculations, we used the overlap Dirac operator. With the 
described setting, we see that, as expected, the scalar density of the zeromode is smeared out in the temporal direction 
for the vortices with bigger temporal extent. The spatial localization, however, is more or less independent of the 
temporal extent of the vortex. Thus, all what really matters is that the vacuum to vacuum transition occurs, not how 
fast it occurs. Let us discuss this in a little bit more detail. First, we note that the scalar density of the zeromode is 
distributed, almost perfectly, spherically symmetric around the midpoint r*o of the vortex. Keeping this in mind, the 
discussion of the localization in 3 dimensional space reduces to a discussion of the localization in \f — r*o|. Thus we 
study 

p{r,t) = j^^2S(r,\r-F \)^o(r,t)M^t) with N(r) = ^ 6(r, \f - f |) . (19) 

r r 

In other words, the quantity p(r, t) stands for the mean value of the scalar density of the zeromode, calculated for 
points belonging to the same \ f — fo\ and t. In Fig. [2] the results for p{r,t) for generalized spherical vortices with 
At = 1 and At = 5 arc compared. One can see the similarity between the result for the two different At. 



VI. DISCUSSION AND OUTLOOK 



In this paper, wc identified the continuum object corresponding to the previously considered spherical vortex as a 
vacuum to vacuum transition in temporal direction. The discrepancy between the gluonic lattice topological charge 
and the index of the lattice Dirac operator, described in previous papers, turned out to be a discretization effect in 
the temporal direction. Starting from the continuum spherical vortex, we constructed a generalized spherical vortex 
on the lattice. We demonstrated the similarity to the original spherical vortex by using fermions as probes. We also 
showed that, with an appropriate choice of parameters, the action of the generalized spherical vortex can be quite 
small, as small as about 5/3 of the one- instanton action. For more details see pH ]. 

It is known that topological charge contributions from center vortices, due to intersection (and writhing) points, 
emerge from vortex structures lying in a single U(l) subgroup (22j . Color rotations in such structures are suppressed by 
the action. For example, intersections of vortices with orthogonal color structures give maximally negative plaquettes 
and are thus suppressed in the continuum limit. However, this does not mean that the whole vortex structure can 
be gauge transformed to a single U(l) subgroup. We conjecture, for entropic reasons, that in the confined phase the 
color structure of center vortices contributes to the topological charge and chiral symmetry breaking. 

Further investigations should therefore aim at quantifying these kinds of contributions in Monte Carlo generated 
configurations. First measurements show that the total vortex surface definitely covers the full S 3 , i.e. vortex 
plaquettes are distributed uniformly among the entire SU(2) color palette. However, color vectors are not gauge 
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FIG. 2. p(r,t) as defined in Eq. (|19|) with normalization p max = 1. The configuration investigated is given by Eq. (| 18p with 
At = 1, (red) diamonds, and At = 5, (blue) circles. The parameters of a(r) are given by R = d = 6 in both cases. The 
calculations have been performed on a 24 3 x 12 lattice. The different curves belong to different values of t. For At = 1 we have 
t e {1 = 2, 3 = 12, . . . , 7 = 8} and for At = 5 we have t G {3 = 4, 5 = 2, . . . , 9 = 10}. The highest curve corresponds to the first 
pair of t-values, the lowest curve to the last pair of t- values. In between the curves fall monotonically. The solid (magenta) line 
represents the average over t for At = 1, the dash-dotted (green) line represents the same for At = 5. 



invariant and therefore the number of full coverings is not so easy to measure, even though it is a gauge invariant, 
topological quantity. Further, writhing points dominate the topological susceptibility (23lI24|. even in 5(7(3) [25j, and 
therefore color contributions might only play a sub-dominant role like intersection points. But the colorful spherical 
vortex may act as an ansatz for a model of chiral symmetry breaking, as it shows properties similar to instantons pfjj . 
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